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Abstract. We calculate the eigenspectrum of random walks on the Eden tree in two 
and three dimensions. From this, we calculate the spectral dimension ds and the walk 
dimension and test the scaling relation ds = 2df/dyj (= 2d/dyj for an Eden tree). 
Finite-size induced crossovers are observed, whereby the system crosses over from a 
short-time regime where this relation is violated (particularly in two dimensions) to 
a long-time regime where the behavior appears to be complicated and dependent on 
dimension even qualitatively. 
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1. Introduction 



Treelike structures arise in many situations in statistical physics, including clusters 
formed by irreversible growth processes modeled, e.g., by diffusion-limited aggregation 
[1] (which are treelike on large scales) and dielectric breakdown [2]. In this paper, we 
consider a particularly simple tree structure called Eden tree [3] which is formed on 
a lattice by a modification of the usual Eden process [4] (in which a compact cluster 
grows) so that an empty site which is a neighbor to more than one occupied site 
becomes ineligible for occupation. 

The study of diffusion on treelike structures is interesting in that it may present 
qualitatively different behavior compared with the structures with many loops 
(particularly those that connect large branches). The Eden tree model was studied 
in detail by Dhar and Ramaswamy [3] using various methods. They calculated the 
spectral dimension ds [5] for two- and three-dimensional Eden trees by applying the 
node counting theorem [6] and also the random walk dimension d^ [7] by a direct 
simulation of random walks as well as through the backbone statistics of the tree, and 
discussed them in terms of a scaling argument. A surprising conclusion was drawn 
from these discussions; they concluded that a typical long walk samples only order-one 
segments of the backbone (including the dangling ends attached to this segment) so 
that the usual scaling relation [8] 

ds = 2df/dy, (1) 

does not hold. (Here df — d since the Eden trees are compact.) While the 
calculation of ds (particularly in two dimensions) seemed accurate, the graphs from 
which d^ was deduced appeared to show considerable curvature. Thus a more accurate 
determination of these exponents appeared desirable to test the violation of relation 
(1) better. 

The scaling law (1) is expected to be valid in general because it follows by assuming 
only: (a) the equivalence of the vibrational and diffusional problems (Alexander 
and Orbach [8]) and (b) that a random walk samples essentially all sites more or 
less uniformly. (In this paper, we shall loosely call the latter condition ergodic for 
simplicity.) In fact this relation is sometimes used as the definition of the spectral 
dimension ds- In their work [3], Dhar and Ramaswamy essentially reject assumption 
(b) thereby rejecting the relation for the Eden tree. 
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This is surprising because it appears at first that any random walk satisfying the 
detailed balance condition should in principle be ergodic (in the way we use the term in 
this paper). This should be true independent of the structure and thus for example (1) 
holds for Euclidean lattices in all dimensions and also for critical percolation clusters 



In this paper, we reexamine this rather curious observation of Ref.[3], using a 
method which allows a more accurate numerical evaluation of both dg and dyj from 
the same data. This is the method of calculating the eigenspectrum of the so-called 
hopping probability matrix W by Saad's diagonalization method [10] and making use 
of the Laplace transform relationship between this spectrum and quantities such as the 
return-to-origin probability and the velocity autocorrelation function [11,12]. (Here 
W has elements Wij which are just the probability for a random walker at site j to 
hop to site i in the next time step.) This corresponds to averaging over all random 
walks starting from every point of the cluster. The method was used previously [12,13] 
for two- and three-dimensional percolation clusters very successfully and we expect it 
to work just as well for the Eden trees. 

Thus we first construct an Eden tree of size S by simulation on the square and 
simple cubic lattices, create the matrix W for a random walk (we use a blind ant 
model where the random walker has an equal probability to hop between any pair 
of neighboring sites), and calculate the eigenvalues and eigenvectors of W for the 
largest M (typically a few hundred) eigenvalues to very high precision (typically up 
to 6 decimal places) using Arnoldi-Saad algorithm [10]. The eigenvalue density n(A) 
is then expected [13] to scale as: 



One can also construct another interesting function 7r(A) = n{X)ax{X — 1)^ [14] from 
the eigenvalues and eigenvectors of W, where ax are the coefiicients entering the mean 
position autocorrelation < r{t) ■ r(0) >= axX^. The coefficients ax are determined 
when the stationary initial distribution is expanded in terms of the eigenvectors of W. 
This function is expected [12,13] to scale as 



[9]. 



n(A) ~ I In A 



^3/2-1 



(2) 



7r(A) ~ I In A 



1-2/d, 



(3) 
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2. Numerical results 



Our numerical results for n(A) and 7r(A) are plotted in Fig.l and Fig. 2 respectively. 
(The values of tt plotted have been scaled up by a factor of the size S of the tree 
compared with the convention used in Ref.[13] in order to separate different S in the 
figure.) The cluster sizes for which the calculations were made are S = 2500, 5000, 
10000 for the square lattice for which 200, 400, and 200 independent realizations were 
averaged over (respectively), and S = 2500, 5000 for the simple cubic lattice for which 
200 and 400 independent realizations were averaged (respectively). 

From Fig. 1(a) for the square lattice, there is apparently a sharp cluster-size 
dependent crossover which divides a region of larger | lnA| (corresponding to shorter 
time) where an excellent power-law fit can be made, from the smaller |lnA| region 
where another power-law fit may also be possible. We will see later that the locations 
of these crossovers are consistent with the finite size scaling with the walk dimension. 
Fig. 1(b) for the simple cubic lattice gives a much less clear indication of such a 
crossover, but nonetheless, the sudden flattening of the data in the regions of smallest 
I In A I is consistent with a crossover interpretation. We observe that the right-hand side 
regions in these figures should give the asymptotic exponent corresponding to Eq.(2) 
as S increases (because this is the region before the finite size effects appear to set in) . 
However, the apparent power-laws in the left-hand side region are still ununderstood. 
Moreover, it is even possible that the qualitative behavior of these crossovers changes 
when S is made much larger; since in the present analysis we are limited to relatively 
small Eden trees, we cannot be sure if there is not another true asymptotic behavior. 

Estimating the slopes of the log-log plots in the right-hand side regions and using 
Eq.(2), we obtain 

4 = 1.22 ±0.02 (d = 2) (4) 

and 

4 = 1.32 ±0.02 (cZ=3), (5) 

where the errors quoted are mainly the fitting errors but also include some error due 
to the choice of the region to fit and that due to the small variation for different cluster 
sizes. These estimates are consistent with those of Ref. [3] but improves the accuracy 
especially in three dimensions. In particular, we rule out the possibility that ds is the 
same for two and three dimensions. 
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The results for 7r(A) in Fig. 2 show an even more interesting cluster-size dependent 
anomaly. In Fig.2(a) for the square lattice, starting from the larger | In A| (i.e., shorter 
time) region, there is a range of about one to two decades where an excellent power 
law fit can be made, and then a sudden and large decrease occurs. This anomaly 
is reproducible in independent batches of clusters and thus not caused by simple 
statistical fluctuations. In Fig. 2(b) for the simple cubic lattice, another sharp, but 
very different type of crossover is observed, where the small |lnA| (or long time) 
region also gives a power-law but with a different exponent. 

Let us first discuss the results shown in Fig.2(a) for the square lattice in more 
detail. Here, if only the flat region for larger |lnA| is fitted to a power-law, the 
exponent comes out to be about 0.29 ±0.01, 0.29 ±0.01, and 0.19 ±0.01 for S = 2500, 
5000, and 10000, respectively. We note that the last value would yield dyj = 2.47±0.03, 
which would be barely consistent with the value dyj = 2.54 ± 0.04 obtained in Ref.[3] 
from the backbone statistics. This also means that these data would violate (1) very 
strongly, since 2d/d^ would be about 1.62 ± 0.02 using this value of dyj. Moreover, 
since the crossover point where the drop begins moves to smaller |lnA| for larger 
cluster size S, we would expect this exponent to be the correct asymptotic exponent 
dyj for S — > oo. 

Our interpretation of this result is as follows: for sufficiently small time (thus 
length) scales (t << r) a typical random walk is trapped (by lack of loops) within 
order-one segment of the backbone (together with its dangling side branches) as 
proposed in [3], which results in a small dyj and the violation of (1). The crossover 
time T should be such that 

T^/'^- oc L (6) 

where L is the length scale of the cluster and d^ is the walk dimension in this regime. 
Translated to the corresponding value in | lnA|, the crossover point should scale as 

llnAil ocT"^ oc L"^-. (7) 

However, as the walk goes beyond this time scale, it is forced by finite size effect to 
sample more and more segments of the backbone. Thus over long times, an effective 
dw increases, leading to the downward tendency in Fig. 2(a) beyond the crossover. Of 
course, for even longer times, the size of the visited region becomes limited by the 
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total size of the cluster and the walk visits every site of the cluster; however, such a 
long-time regime is not reached in our figure. 

Numerically, the crossover values |lnAi|for n{X) are approximately 5 x 10~^, 
2 X 10~^, and 1 x 10"'* for S = 2500, 5000, and 10000, respectively. In comparison, 
L"^- with L = y/S/2 and the estimate d^j = 2.47 gives 3.5 x 10"^, 1.5 x 10"^, and 
6.4 X 10~^, in fair agreement with the discussion above. The discrepancy of about a 
factor of 1.5 in the absolute numbers should be due to the fact that our method is 
equivalent to summing over all random walks starting from all sites of the cluster and 
not just from the seed site at the center. 

Having identified dw, we can return to the crossovers observed in Fig. 1(a) for n(A) 
in two dimensions. These crossovers occur roughly around | lnA| = 1 x 10~^, 5 x 10~^, 
and 2 x 10~^ for 5" = 2500, 5000, and 10000, respectively. The ratios of these numbers 
are again in fair agreement with those from Eq.(7). 

The situation for 7r(A) in three dimensions is much less clear. In Fig. 2(b), 7r(A) is 
plotted for S = 2500 and 5000 for the simple cubic lattice. Similarly to the square 
lattice, there is a flat region for larger | In A| which can be fitted to a power-law and 
what seems to be a crossover to another power-law at smaller | In A|. The first region 
gives slopes of 0.49 ± 0.01 and 0.51 ± 0.01 for S = 2500 and 5000, respectively. This 
would translate to d^ of about 3.92 ±0.08 and 4.08 ±0.09, respectively, and thus, with 
df = d = 3, we would have 2df /d^ of around 1.5 which would yield a much smaller 
deviation from Eq.(l) than in two dimensions. 

Beyond the crossover in the smaller | lnA| region, we have another power-law with 
an exponent of about 0.32 and 0.29 for S = 2500 and 5000, respectively. Since these 
exponents would lead to smaller values of d^, in this region, we would have an even 
stronger violation of the scaling relation (1). 

The crossover point is at about | In Ai| = 1 x 10~^ and 4 x 10~^ for S = 2500 and 
5000, respectively. The ratio of these two numbers is in good agreement with Eq.(7) 
when a value of dw = 4.08 is used. Thus, we tentatively conclude that the asymptotic 
behavior for an infinite Eden tree in three dimensions may also be nonergodic where 
Eq.(l) is violated (but less strongly than in two dimensions). However, we do not 
understand the regime for t >> r. One speculation would be that, in the latter 
regime, the random walk does not become ergodic over the whole cluster but becomes 
trapped in some parts of it, e.g., in the surface of the Eden tree. The walk dimension 
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of this trapping region would then determine the slope after the crossover. However, 
the surface of the Eden model is a particularly diflBcult problem from the numerical 
point of view, with a very slow and non-monotonic convergence to the asymptotic 
behavior [15]. It is also possible that the cluster sizes considered here are not yet in 
the asymptotic regime and that for much larger clusters a different picture emerges. 
Unfortunately no suitable technique is available to calculate ds and accurately on 
very large clusters. 

3. Conclusion and discussion 

In this work wc have calculated the eigenspectrum of the random walk hopping 
probability matrix W and from this computed the critical exponents ds and du, for the 
Eden tree in two and three dimensions. By more accurate calculations, we confirm 
the observation of Ref.[3] that the random walk is not ergodic on two dimensional 
Eden tree and moreover our data suggest that the same is true for three dimensions 
(although less strongly). Here, we use the term ergodic to express the notion that a 
typical random walk visits essentially all sites of the cluster (i.e., at least sites that 
have the same fractal dimension as the entire cluster) in the long time limit. 

A natural question is why the detailed balance in this problem docs not guarantee 
ergodicity. This is understood when one makes a distinction between ergodicity as we 
used the term here and a different statement that any initial distribution eventually 
evolves into the stationary one. The latter statement is true when detailed balance is 
satisfied, but this is a statement about the distribution and not about an individual 
random walk. Thus it is not necessary for an individual random walk to sample 
essentially all sites of the cluster with probability one even in the t ^ oo limit. Rather 
it can so happen that a typical walk is very anisotropic, sampling only order-one 
segment of the tree, but when we consider all possible walks, the entire cluster is 
sampled, if one waits long enough. 

It is also interesting to consider whether a similar situation arises in tree structures 
in general. A typical treelike structure encountered in statistical physics is the 
diffusion limited aggregate (DLA) [1]. The available numerical results [16] from direct 
simulation of random walk displacements indicate that d^i = 2. 56 ±0.10 in d = 2 and 
3.33 ± 0.25 in d = 3, while the simulation of the return-to-origin probability in the 
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same reference gives dg = 1.20 ± 0.1 in d = 2 and 1.30 ± 0.1 in d = 3. Using the 
known estimates [5] of df = 1.68 ± 0.05 in cZ = 2 and 2.5±0.06in(i = 3 together with 
the estimates for from [16], we obtain the value of 2df/dy, to be about 1.31 ±0.1 
in d = 2 and 1.50±0.1 in d — 3. While these values are not inconsistent with the 
direct estimates of ds from [16], they are sufficiently different (particularly in (i = 3) 
to warrant further investigation. 

Indeed, it may be possible that more or less all treelike structures force a typical 
random walk to be anisotropic. If this were the case, it would be of great interest to 
develop a systematic theory of this anisotropy, perhaps with different coherence length 
exponents characterizing the radial and tangential directions in a manner somewhat 
similar to directed percolation [9]. This is, however, clearly out of scope for this work 
and must await further research. 
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Figure captions 



Figure 1. (a) The eigenvalue density n(A) for the Eden tree on the square lattice. 
The data points Q), x , and A correspond to the average over 200, 400, and 200 clusters 
of size S = 2500, 5000 and 10000, respectively. The solid lines are obtained by least 
squares fitting to the data for | lnA| larger than its crossover value (number of points 
used in the fits being 14, 19, and 19 for the three sizes respectively), and the dotted 
lines are from the similar fits to the data smaller than the crossover point; (b) n{X) 
for the simple cubic lattice where the symbols Q ^tnd x correspond also to = 2500 
and 5000, respectively. 

Figure 2. (a) The function 7r(A) for the Eden tree on the square lattice. The symbols 
O, x, and A have the same meaning as in Fig.l; (b) Same for the simple cubic lattice. 
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